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Abstract. The observed average lifetime of the population of star clusters in the Solar Neighbourhood, the Small 
Magellanic Cloud and in selected regions of M51 and M33 is compared with simple theoretical predictions and 
with the results of TV-body simulations. The empirically derived lifetimes (or disruption times) of star clusters 
depend on their initial mass as tdi a cmp oc M c i 0,60 in all four galaxies. TV-body simulations have shown that the 
predicted disruption time of clusters in a tidal field scales as idis pred oc t^ 5 t^ 25 , where t r h is the initial half-mass 
relaxation time and t CT is the crossing time for a cluster in equilibrium. We show that this can be approximated 
accurately by £dis prcd oc M^i 62 for clusters in the mass range of about 10 3 to 10 6 Mq, in excellent agreement with 
the observations. Observations of clusters in different extragalactic environments show that tdis also depends on 
the ambient density in the galaxies where the clusters reside. Linear analysis predicts that the disruption time 

1/2 

will depend on the ambient density of the cluster environment as tdis oc p am (, . This relation is consistent with 
TV-body simulations. The empirically derived disruption times of clusters in the Solar Neighbourhood, in the SMC 
and in M33 agree with these predictions. The best fitting expression for the disruption time is 

t dis = C cnv (A/ cl /lO 4 Af ) a62 (p amb /M pc- 3 )- a5 

where M c i is the initial mass of the cluster and C cnv ~ 300 — 800 Myr. The disruption times of star clusters in M51 
within 1-5 kpc from the nucleus, is shorter than predicted by about an order of magnitude. This discrepancy 
might be due to the strong tidal field variations in M51, caused by the strong density contrast between the spiral 
arms and interarm regions, or to the disruptive forces from giant molecular clouds. 

Key words. Galaxy: open clusters - Galaxy: solar neighborhood - Galaxies: individual: M33 - Galaxies: individual: 
M51 - Galaxies: individual: SMC - Galaxies: star clusters - 



1. Introduction 

The age distribution of star clusters in the disk of the 
Milky Way, can only be explained if galactic clusters dis- 
rupt on a time scale on the order of a few times 10 8 years 
(Oort 1957; Wielen 1971, 1988). Clusters in the LMC and 
SMC survive longer than those in the solar neighborhood 
(e.g. Elson & Fall 1985,1988; Hodge 1987). 

The disruption time of clusters is expected to depend 
on both the internal cluster conditions, such as the initial 
mass, density and velocity dispersion and the stellar initial 
mass function (IMF), and on the external conditions, such 
as the orbit in the galaxy and tidal heating by encoun- 
ters with for example giant molecular clouds. Recently, 
Boutloukos & Lamers (2003, hereafter called BL03) have 
derived an empirical expression for the disruption time 
of clusters as a function of their initial mass in selected 
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regions in four galaxies. They found large differences in 
the disruption times between these regions. This allows a 
crucial test of the theoretical disruption times of clusters 
predicted with TV-body simulations. 

The purpose of this paper is three-fold: 

(a) to explain the empirical dependence of the cluster dis- 
ruption times on cluster mass; 

(b) to explain the strong dependence, found by BL03 of 
the disruption time on the conditions in the host galaxies; 

(c) to confront the predicted disruption theory with the 
empirically derived disruption relations. 

In Sect. 2 we describe the predicted dependence of the 
disruption times of star clusters on their initial mass and 
on their environment. In Sect. 3 we review the method 
used for the determination of the disruption times based 
on a statistical analysis of large cluster samples and the re- 
sulting expressions for the disruption times in four galax- 
ies. In Sect. 4 we discuss the observed dependence of the 
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disruption times on the environment of the clusters in their 
host galaxies and we compare it with the predictions. We 
discuss the peculiar case of the clusters in M51. The con- 
clusions are in Sect 5. 

2. The predicted dependence of the disruption 
time on cluster parameters and environment 

The evolution of an isolated star cluster is characterized by 
three distinct phases: (1) infant mortality (the dissolution 
of unbound clusters) within about 10 Myr (Bastian et al. 
2004; Whitmore 2004), (2) stellar evolution dominates the 
cluster mass loss in the first ~ 10 s years (see Takahashi & 
Portcgics Zwart 2000) followed by (3) a relaxation domi- 
nated phase in which the cluster mass loss is driven by its 
internal dynamical evolution and the external influence of 
the tidal field of the host galaxy. If the cluster survives the 
first two phases, its lifetime will be dominated by the third 
phase as it generally lasts much longer, until complete dis- 
solution. During this episode mass loss is about constant 
with time and appears to depend weakly on the density 
profile of the cluster (Spitzer 1987; Portegies Zwart et al. 
2001b, but see Portegies Zwart et al. 1998, for some com- 
plications near the end of the lifetime of the star cluster). 
During the third phase there are two important effects 
which drive the cluster dissolution. These are internal two- 
body relaxation and the interaction between the cluster 
and the tidal field of the host Galaxy. 

Internal two-body relaxation drives the cluster evapo- 
ration on a very long time scale of many initial half- mass 
relaxation times (Baumgardt et al. 2002). The half- mass 
relaxation time scale can be written as (Spitzer 1958): 
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where N is the initial number of stars in the clusters, 
(m) = M c i/N is the mean mass of the stars and r^m is 
the half-mass radius. The Coulomb logarithm A is ap- 
proximately proportional to N (see Gierz & Heggie 1996; 
Spinnato et al. 2003). 

The mass loss rate for clusters in a tidal field is much 
higher and the lifetime shorter than for clusters in isola- 
tion. Baumgardt (2001) found that for clusters in a tidal 
field the disruption time depends not only on the relax- 
ation time but also on the crossing time t cr because, due 
to their circulation within the cluster, stars with sufficient 
energy still need a considerable time to reach the outskirts 
of the cluster from where they can escape. He showed that 
the disruption time can be approximated as 
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where the crossing time for a cluster in equilibrium is 
icr « r 3 h ^/y/GM^. (3) 
and so t cr oc i r h In A/N. 



Baumgardt (2001) showed on theoretical arguments 
that x is expected to be about 3/4. A more accurate esti- 
mate of the value of < x < 1 can be found empirically 
from A^-body simulations (see below). Combining Eqs.^[21 
andean expression for the disruption time can be derived 
of the form 
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where we replaced A with 7 AT, with 0.01 < 7 < 0.4. 

Baumgardt & Makino (2003, hereafter BM03) have 
performed an extensive set of simulations of clusters in 
the Galactic halo tidal field and fitted the results to find 
the constant C and the value of x in the approximation 
of Eq. They used a logarithmic potential of the form 
0(i?c) = VqIii(Rq), where Rq is the distance to the 
galactic center and Vq is the circular velocity. Their clus- 
ters started with a tidal radius equal to the tidal radius 
of the external field, 
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If we assume the tidal radius scales linearly with the half- 
mass radius, which holds for clusters with the same density 
profile, Eq. [SJcan be combined with Eq. which yields 
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where (3 is a constant whose value can be found empiri- 
cally. 

BM03 fitted this relation to the results of their simu- 
lations of clusters at different distances from the Galactic 
center. For clusters that initially followed a King (1966) 
profile with central concentration Wo = 5.0 they found (i 
= 1.91 and x = 0.75. For clusters with W = 7.0 they 
found (3 = 1.03 and x = 0.82. For 7 a value of 0.02 was 
adopted. We found that the factor (3{N/ ln(0.02N)} x in 
Eq. [5] for both combinations of /3 and x can be approxi- 
mated by a power law 1.95A' ' 62 with high accuracy. This 
is illustrated in Fig. Q 

The density of the clusters depends on the ambient 
density which depends on the adopted potential field. The 
ambient density, for the potential field adopted in the 
study by BM03 can be found by applying Poisson's law 
which gives 
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Combining this with Eq. and assuming an initial mean 
stellar mass of (m) = 0.54M Q , as done by BM03 based 
on the stellar IMF of Kroupa (2001), we can write the 
disruption time as a function of cluster mass and local 
ambient density: 
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Fig. 1. The function f3{N/ \n(0.02N)} x of equation for 
the two combinations of f3 and x for clusters with Wo = 5 
(solid line) and 7 (dashed) from BM03 and the power law 
approximation 1.95 N om7 (dots). 

with M c \ in units of Mq and p a mb in Mqpc -3 . The con- 
stant indicates the time when 95 per cent of the initial 
cluster mass is lost. If the initial mean stellar mass is 0.70 
M Q , predicted by the stellar IMF from Scalo (1986), the 
constant decreases from 810 to 690. (For clusters in el- 
liptical orbits with ellipticity e moving within a logarith- 
mic potential field the disruption time is shorter by factor 
1 — e, if Pamb is the density at the apogalactic radius.) The 
positive exponent 0.62 indicates that the disruption time 
increases with mass, as expected. 

The results of TV-body simulations and their predicted 
dependence on the ambient density are plotted in Fig. [21 
These are based on iV-body simulations of clusters of 10 4 
Mq from different authors. These simulations are from 
Portegies Zwart et al. (1998; 2002), and BM03 . All simula- 
tions were performed with NB0DY4 (Aarseth 1999; BM03) 
and the Starlab (Portegies Zwart et al. 2001a) 1 software 
environment with both codes running on the GRAPE fam- 
ily of special purpose computers (Makino 2003). The mod- 
els by Portegies Zwart et al. (2002) were calculated for star 
clusters near the Galactic center, i.e. they are representa- 
tive for a high density ambient medium. 

The dashed line shown in this figure represents the 
mean relation derived by BM03 and described by Eq. (jSJ) . 
This line passes through the data of their individual mod- 
els. 

The full line is the mean relation predicted by Portegies 
Zwart et al. (2001b). It is derived from the argument that 
the lifetime of a cluster will depend on the two-body re- 
laxation time near the tidal radius as idis — 0.3i r t- The 
tidal radius can be obtained from Eq. [SJ but the required 
average cluster density p c \ is not readily available from the 
observations. Since the cluster is in equilibrium with its en- 
vironment one can estimate p c i/pamb — 3 for a standard 
Roche-solution by assuming that the clusters follow a King 
(1966) model and that the tidal radius of the King model 

1 see also http: //www .manybody. org 



equals the Jacobi radius of the star cluster in the tidal 
field of the Galaxy (see e.g. Binney & Tremaine 1994). 
This predicts that the constant in Eq. ((SJ) is about 330 
Myr. This is 0.4 times as small as the value derived from 
the models by BM03 . The line passes through the data 
of the individual models of Portegies Zwart et al. (1998, 
2002). 

The simulations by Portegies Zwart et al. (1998, 2002) 
result in shorter disruption times than the simulations by 
BM03, by about a factor of 2.5. This difference can be ex- 
plained by variations in the initial conditions and between 
the two codes. The initial mass function used in Portegies 
Zwart et al. (2002) is taken from Scalo (1986) and has a 
mean stellar mass of (to) ~ 0.7Af Q . BM03 adopted the 
initial mass function of Kroupa (2001), which has a mean 
stellar mass of 0.54 Mq. Portegies Zwart et al. (1998) and 
BM03 use King (1966) models for the initial density dis- 
tribution, whereas Portegies Zwart et al (2002) adopt self 
consistent Heggie & Ramamani (1995) models. But also 
the stellar evolution in the calculations are slightly dif- 
ferent: the stellar evolution in BM03 is based on Pols et 
al. (1998) whereas the calculations of Portegies Zwart et 
al. use the stellar evolution predicted by Eggleton et al. 
(1989). Moreover, the end of the cluster lifetime is defined 
differently in the various calculations: BM03 and Portegies 
Zwart et al. (2002) uses the time when 5 percent of the 
initial mass remains, whereas Portegies Zwart et al. (1998) 
use the time when 10 percent of the initial mass remains. 
The most important difference, however, is in the treat- 
ments of the tidal field and the way stars escape from the 
cluster. Portegies Zwart et al. (1998) did not include a 
full tidal potential, but a simple cut-off radius. This easily 
leads to a shorter disruption time for the cluster. BM03 
and Portegies Zwart et al (2002), on the other hand adopt 
a self-consistent tidal field up to the quadrupole moment. 2 

A few notes about these estimate: 

(a) The constant of expression @ was derived for clusters 
moving in circular orbits in the Galactic potential field. 
Clusters in elliptical or inclined orbits will experience 
disk or bulge shocking that reduces their lifetime (Spitzer 
1987, Zhang & Fall 1999; BM03). However the propor- 
tionality with mass and density remains the same. 

(b) Encounters with giant molecular clouds (GMCs) or 
spiral density waves are not included in these calculations. 
In situations where these effects are likely to occur, e.g. 
for clusters in orbits in a galactic plane, the constant will 
be smaller than estimated above. 

Based on these considerations we can expect the fol- 
lowing expression for the disruption times of star clusters 



2 In an extensive collaborative experiment Heggie et al. 
(1998) perform a detailed comparison between the N-body 
codes used by Portegies Zwart and Baumgardt (and others). 
The difference in cluster lifetime is between 20 and 40%. (see 
http: //www . maths . ed.ac.uk/~douglas/ experiment .html ) 
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Fig. 2. The predicted dependence of the disruption time, 
ti, of clusters of initial mass of 10 4 M on the mean am- 
bient density, Pamb- The symbols refer to iV-body simula- 
tions; filled squares are from Portegies Zwart et al. (1998, 
from left to right models IR16, IC16 and FH16), filled tri- 
angles are calculations for star clusters near the Galactic 
center from Portegies Zwart et al (2002, from left to right 
models R150W4, R90W4 and R34W4). The open circles 
represent models from BM03 (from left to right models 
V, I, II and IV with 16384 stars). The dashed line is the 
relation predicted by BM03, Eq.|S| and the full line is the 
relation predicted by Portegies Zwart et al. 

due to tidal interactions in different extragalactic environ- 
ments 

tdis = Ccnv [wif^ ) (9) 

where C cnv ~ 300 — 800 Myr. If encounters with giant 
molecular clouds or disk shocking (for clusters in orbits 
tilted with respect to the galactic plane) becomes impor- 
tant, the disruption times might be shorter than given by 
Eq. @. 

3. The empirically determined disruption times 

The disruption times of systems of star clusters in dif- 
ferent extragalactic environments have been determined 
statistically by BL03 from samples of star clusters in four 
galaxies. We will summarize their method and results. 

Suppose that the disruption times of clusters in some 
region of a galaxy can be written as 

t dis = n x (Afci/io 4 ) 7 (10) 

where M c \ is the initial cluster mass (in M@) and £4 is the 
disruption time (in yrs) of a cluster with an initial mass 
of M c \ — 10 4 Mq. Assume also that the cluster forma- 
tion rate is constant over the period of time for which ti 
and 7 are determined, and that all clusters had the same 



initial stellar IMF. Then, the values of 7 and ti can be 
derived from the mass and age histograms of cluster sam- 
ples with a well determined brightness limit. Boutloukos & 
Lamers (BL03) have applied this method for an empirical 
determination of 7 and ti in specific regions of M51, M33, 
(observed with HST-WFPC2) and in the SMC and the 
solar neighbourhood. The method has been improved by 
Lamers (2004) and by Gieles et al. (2004c), who showed 
that the distribution of the observed clusters in a mass 
versus-age-diagram provides a more powerful method for 
deriving ti and 7 than the mass and age histograms sep- 
arately. 

The results are given in Table 1. Column 2 gives the 
range of galactocentric distances of the cluster sample. 
Columns 3, 4 and 5 give the number of clusters and their 
age and mass range. Columns 6 and 7 give the resulting 
values of the disruption time ti and 7. Column 8 gives 
the ambient density of the host galaxy at the location 
of the cluster samples used for the determination of the 
disruption time (see Sect. 4). 

In the determination of 7 we adopted a cluster IMF 
with a = 2.0 (Zhang & Fall 1999, Bik et al. 2002, de 
Grijs et al. 2003). The values of t 4 and 7 of the SMC is 
from BL03. The values for M51 have been redetermined 
by Gieles et al. (2004c) on the basis of the extended cluster 
sample of Bastian et al. (2004) and the improved method 
mentioned above. The values for M33 have been redeter- 
mined by Lamers et al. (2004b), on the basis of the ex- 
tended cluster sample by Chandar et al. (2001, 2002). The 
new value of ti for this galaxy is considerably higher than 
originally derived by BL03, due to the largely improved 
cluster sample and a redetermination of the cluster ages. 
The values for the Solar neighbourhood have been rede- 
termined by Lamers et al. (2004b), based on the clusters 
sample by Loktin et al. (1994). 

Notice that the values of 7 of the four galaxies are the 
same, within the uncertainty, but the values of ti differ 
strongly. The mean value of the empirically determined 
value of 7 and its uncertainty is 

<7om P ) = 0.60 ±0.02. (11) 

where 0.02 is the uncertainty in the mean value, which is 
much smaller than the uncertainty in the individual mea- 
surements. This exponent agrees very well with 7 pre d = 
0.62 predicted by BM03, Eq. [8] on the basis of their nu- 
merical simulations. The agreement between the observed 
and the predicted value of 7 was first pointed out by Gieles 
et al. (2004d), but based on slightly different theoretical 
arguments. 

4. The dependence of the disruption time on the 
conditions in the host galaxy 

Table 1 shows that the constant ti, which is the disruption 
time of a clusters of M c \ — 10 4 M Q , differs greatly between 
the regions in the galaxies that were studied. This indi- 
cates that the disruption time depends strongly on the 
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Table 1. The parameters of the disruption time: idis = £4 x (Af c i/10 4 ) 7 



Galaxy 


rGal 


Nr 


Age 


range 


Mass 


range 


log £4 


7 


log Pamb 


Remark 




kpc 


clusters 


log 


(yrs) 


log 


M@ 


log (yrs) a 




M Q pc" 3 




M51 


1.0 - 5.0 


1152 


6.3 


- 9.0 


2.6 


- 5.6 


7.85 ±0.22 (7.6) 


0.57 ±0.10 


-0.70 ± 0.30 


2 


M33 


0.8 - 5.0 


147 


6.5 


- 10.0 


3.6 


- 5.6 


8.80 ±0.20 (8.1) 


0.60 ±0.15 


-0.66 ±0.30 


3 


Galaxy (d < lkpc) 


7.5 - 9.5 


184 


7.2 


- 9.5 






8.75 ±0.20 (9.0) 


0.60 ±0.12 


-1.00 ±0.04 


1,4 


SMC 


- 4 


314 


7.6 


- 10.0 






9.90 ±0.20 (9.9) 


0.61 ±0.08 


-2.10 ±0.30 


1 
















Mean 


0.60 ±0.02 
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(a) The value in parenthesis is the original value from BL03 

(1) t± and 7 from age distribution only 

(2) Improved determination of tt and 7 by Gieles et al. (2004c), based on new cluster sample of Bastian et al. (2004) 

(3) Improved determination of £4 and 7 by Lamers et al. (2004b) based on extended cluster samples by Chandar et al.(2001, 
2002) 

(4) Improved determination by Lamers et al. (2004a) based on cluster sample of Loktin et al. (1994) 

(5) The quoted error is the uncertainty in the mean value 



local conditions in the host galaxy. In Sect. (2) we have 
shown that £4 is expected to scale with the inverse square- 
root of the mean density in the host galaxy, if evolutionary 
mass loss, relaxation driven evaporation and tidal strip- 
ping are the dominant mass loss mechanism for clusters. 
In this section we compare this predicted dependence with 
the empirically derived values of the disruption time £4. 

4.1. The ambient density in the host galaxies 

We estimate the ambient density p am b at the galacto- 
centric distances of the observed environments for which 
BL03, Gieles et al. (2004c) and Lamers et al. (2004b) de- 
termined the value of £4. 
M51: 

The clusters are at galactocentric distances of 1 to 5 kpc. 
The column density of M51 decreases exponentially as 
a = 540exp(-r/4.65) Mq pc -2 , with r in kpc (Salo & 
Laurakainen 2000) . The total disk mass within a ring of 1 
to 5 kpc is 2.0 x 10 10 Mq and the mean surface density in 
that annulus is 2.6 x 10 2 Mq pc" 2 . Van der Kruit (2002) 
has shown that spiral galaxies have a vertical scale-height 
of 1/7 times the radial scale- length. This implies a value 
of h z = 0.66 kpc and a mean mid-plane density between 
1 and 5 kpc of 0.20 Mq pc -3 . The structural parameters 
of M51 derived by Athanassoula et al. (1987) result in a 
very similar estimate of the mid-plane density. 
M33: 

The sample of clusters is located at galactocentric dis- 
tances between 1 and 5 kpc. More than 90 percent of the 
clusters are between 1 and 4 kpc. The column density 
of M33 decreases exponentially as a = 525 exp(— r/1.45) 
Mq pc -2 , with r in kpc (Athanassoula et al. 1987). The 
total disk mass within a ring of 1 to 4 kpc is 4.2 x 10 9 Mq 
and the mean surface density is 90 Mq pc -2 . Adopting 
an effective thickness of the disk of 1/7 of the radial scale- 
length, i.e. 0.21 kpc, we find a mean mid-plane density of 
0.22 Mq pc' 3 . 



Solar Neighborhood: 

The mean gas and stellar density in the solar neighbor- 
hood, derived from Hipparcos data is 0.10±0.01 Mq pc -3 
(Holmberg & Flynn 2000). 
SMC: 

The sample of clusters is distributed throughout the SMC. 
With a total mass of about 2x 10 9 Mq (Mathewson & Ford 
1984) and a radius of the central core of 4 kpc (Caldwell & 
Laney 1991) we hnd a mean density of 8 x 10 -3 Mq pc~ 3 . 

The ambient densities are listed in Table 2] where we 
have adopted an uncertainty of a factor 2 in the densities, 
except for the solar neighbourhood. The relation between 
£4 and the ambient densities is shown in Figure |3| The 
figure shows a clear trend of decreasing disruption time 
with increasing ambient density. 

4.2. Comparison between observed and predicted 
disruption time scales 

We compare the observed values of £4 for the different 
galaxies with the predictions based on A^-body simula- 
tions, discussed in Sect. 2. These simulations were done 
for clusters in Galactic environments. However, by using 
the expected scaling laws (Sect. 2), we can apply the re- 
sults of these simulations to the conditions in the galaxies 
where the cluster disruption times were measured. 

In Fig. 0] we compare the observed relation between £4 
and Pamb with the predicted mean relations derived from 
Af-body simulations of clusters of 10 4 Mq from Portegies 
Zwart (1998, 2002) and BM03, shown in Fig.H 

The empirically derived values of £4 roughly follow the 
relations predicted by the A^-body simulations within the 
uncertainty of the data, except for the clusters in M51, 
whose disruption time is almost an order of magnitude 
shorter than expected. 

There are several possible explanations for the short 
cluster lifetimes in M51, some of which originate from 
variations in the birth conditions of the clusters or from 
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Fig. 3. The relation between the measured values of the 
disruption time £4 of clusters with an initial mass of f 4 
Mq and the mean ambient density p a mb in Mq pc~ 3 , 
of the galaxies at the galactocentric distances where the 
studied clusters were located. 
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Fig. 4. Comparison between the observed (large dots) and 
predicted (lines) disruption time £4 of clusters with an ini- 
tial mass of 10 4 Mq as a function of the mean density p a mb 
in Mq pc~ 3 of the host galaxy. The solid curve present the 
predictions by Portegies Zwart et al. (1998, 2002) and the 
dashed line by BM03 . 

star formation histories of the host galaxy. For example, 
the constant used in Eq. (JjJJ depends via (to) , quite sen- 
sitively on the adopted initial mass function. If we adopt 
a Salpeter (1955) initial mass function with a lower mass 
limit of IMq, instead of 0.1 Mq, the mean stellar mass 
increases by about an order of magnitude. Such dramatic 
change to the initial mass function results via Eqs. in 
a reduction of the cluster lifetime by about a factor 3. 
Such a lower mass cut-off for the initial mass function is 



controversial, but is suggested in young star clusters like 
MGG-11 (McCrady et al. 2003) and MGG-F, both in M82 
(Smith and Gallagher 2001) and in several of the young 
star clusters in the Antennae galaxies (Mengel et al. 2002). 

Since M51 experienced a recent star formation event, 
about 60 Myr ago (Bastian et al. 2004), possibly part of 
the discrepancy can be explained by variations in the IMF. 
(This starburst event was taken into account in the deter- 
mination of the disruption time for M51 by Gieles et al. 
2004b). Young star formation environments also tend to 
have higher metalicity, which drives stronger winds and 
therefore more stellar mass loss. Another possible expla- 
nation for the short disruption time of clusters in M51 is 
related to the strength of the tidal field variations. M51 
has the highest density contrast between the arm and in- 
terarm regions of any spiral galaxy (Rix & Rieke 1993), 
so the tidal field variations experienced by star clusters 
will be higher than normal, which might explain the short 
disruption time of the star clusters in M51. This last effect 
is studied by Gieles et al. (2004a). 

5. Conclusions 

We analyzed the disruption time of star clusters in the 
Solar Neighbourhood, the SMC and in selected regions of 
M33 and M51. We find the following result. 

(a) Within each region, there is a clear relation between 
the empirically derived cluster disruption time and the 
cluster mass of the type £dis oc M 1 . The mean empirical 
value of 7emp = 0.60 ± 0.02 agrees very well with the pre- 
dicted value of 7 pre d = 0.62, derived from direct iV-body 
simulations of galactic star clusters. 

(b) Comparing the disruption times in different galaxies, 
we find that the disruption time depends on the ambient 
density in the environment where the clusters reside. N- 
body simulations of clusters at different galactocentric dis- 
tances predict that the disruption time decreases with the 
ambient density as tdis oc /9 a mb -0 ' 5 - The empirical values 
of three galaxies, SMC, M33 and the Solar neighbourhood 
agree roughly with this prediction. However, the disrup- 
tion time of star clusters in the inner 4 kpc of M51 is 
shorter than predicted by about an order of magnitude. 
The difference is most likely due to large tidal field varia- 
tions, due to the high density contrast in the spiral arms 
and interarm regions of M51, or due the starburst trig- 
gered by the last encounter with its companion NGC 5195. 

The disruption time of star clusters of initial mass M c \ 
in circular orbits in an environment with a mean ambient 
density p am b is given by Eq. (JHJ ■ The disruption time may 
be shorter in starburst regions. 
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